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Using large scale numerical simulations we analyze the statistical properties of fracture in the 
two dimensional random spring model and compare it with its scalar counterpart: the random fuse 
model. We first consider the process of crack localization measuring the evolution of damage as the 
external load is raised. We find that, as in the fuse model, damage is initially uniform and localizes 
at peak load. Scaling laws for the damage density, fracture strength and avalanche distributions 
follow with slight variations the behavior observed in the random fuse model. We thus conclude 
that scalar models provide a faithful representation of the fracture properties of disordered systems. 



I. INTRODUCTION 

The statistical properties of fracture in disordered me- 
dia represent an intriguing theoretical problem with im- 
portant practical applications lj . The presence of disor- 
der naturally leads to statistical distributions of failure 
stresses, accumulated damage, acoustic activity, crack 
shapes and so on. The application of a standard contin- 
uum descriprion based on elastic equation cannot capture 
the effect of fluctuations and hence the effect of disorder 
has to be considered explicitly. A well established way to 
deal with this problem relies on lattice models, in which 
the medium is described by a discrete set of elastic bonds 
with randomly distributed failure thresholds 0. In the 
simplest approximation of a scalar displacement, one re- 
covers the random fuse model (RFM) where a lattice of 
fuses with random thresholds are subject to an increasing 
external current 0- 

The RFM has been extensively investigated in the 
last twenty years, mainly using numerical simulations 
[B 111 III El El El 01 Q The type of behavior at macroscopic 
fracture is significantly influenced by the amount of dis- 
order Q. When the disorder is narrowly distributed, ma- 
terials breakdown without significant precursors. As the 
disorder increases, substantial damage is accumulated 
prior to failure and the dynamics resembles percolation 
|9|. Indeed, in the limit of infinite disorder, the dam- 
age accumulation process can exactly be mapped onto 
a percolation problem E3- It has been suggested that 
for strong, but finite, disorder fracture should be inter- 
preted as a first order transition near a spinodal point 
, 6| . In addition, the fracture of the RFM is preceded by 
avalanches of failure events 0, 0> E3- These are 
reminiscent of the acoustic emission activity observed in 
experiments and their distribution follows a power law. 
Finally, the RFM has also been used to compute the 
fracture str eng th distribution and the related size effects 
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Modeling the elastic medium using the RFM intro- 
duces drastic approximations in terms of the discretiza- 
tion process, quasistatic dynamics and the scalar nature 
of the interactions. It is thus important to clarify if the 
observations made in the RFM carry over to more com- 
plex and realistic situations. In this paper, we address the 
problem of the scalar (electric) interactions of the RFM, 
by comparing it with a tensorial central force model, the 
random spring model (RSM) E3 ■ The model is a tenso- 
rial counterpart of the RFM: it has quasistatic dynamics, 
random thresholds, but fuses and currents are replaced 
by elastic springs and forces. Dynamic effects have been 
instead considered in Refs. |20tl21j. 

After discussing the model in Sec. II, we consider the 
typical statistical measures performed using the RSM: 
damage localization and average damage profiles are re- 
ported in Sec. Ill, while mean damage scaling and dam- 
age distributions are discussed in Sec. IV. and Sec. V 
respectively. In Sec. VI and VII we discuss the frac- 
ture strenght distribution and the size effect on the mean 
strength. The avalanche behavior is analyzed in Sec. VII 
and a summary is reported in Sec VIII. We have not 
analyzed the roughness of the final crack since in sev- 
eral instances the spring networks fail because of loss of 
rigidity. 



II. THE RANDOM SPRING MODEL 

In the RSM, the lattice is initially fully intact with 
bonds having the same stiffness, but the bond breaking 
thresholds, t, are randomly distributed based on a thresh- 
olds probability distribution, p(t). The bond breaks ir- 
reversibly, whenever the force in the spring exceeds the 
breaking threshold force value, t, of the spring. Peri- 
odic boundary conditions are imposed in the horizontal 
direction and a constant unit displacement difference is 
applied between the top and the bottom of lattice system. 
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Numerically, a unit displacement, A = 1, is applied at 
the top of the lattice system and the equilibrium equa- 
tions are solved to determine the force in each of the 
springs. Subsequently, for each bond j, the ratio between 
the force fj and the breaking threshold tj is evaluated, 
and the bond j c having the largest value, maxj-^, is ir- 
reversibly removed. The forces are redistributed instan- 
taneously after a bond is broken implying that the stress 
relaxation in the lattice system is much faster than the 
breaking of a bond. Each time a bond is broken, it is nec- 
essary to re-equilibrate the lattice system in order to de- 
termine the subsequent breaking of a bond. The process 
of breaking of a bond, one at a time, is repeated until the 
lattice system falls apart. For the RSM, we consider a tri- 
angular lattice system network and a uniform probability 
distribution for thresholds disorder, which is constant be- 
tween and 1. The diamond lattice (square lattice with 
bonds inclined at 45 degrees) spring system exhibits cer- 
tain unstable modes and hence is not considered. Figure 
n presents the envelope of a typical force-displacement 
response obtained using the RSM. The peak load of the 
lattice system is defined as the maximum force of the 
force-displacement response. 



Numerical simulation of fracture using large lattice 
networks is often hampered due to the high computa- 
tional cost associated with solving a new large set of 
linear equations every time a new lattice bond is bro- 
ken. In this study, we use the multiple-rank sparse 
Cholesky factorization updating algorithm developed in 
Ref. for simulating fracture using discrete lattice sys- 
tems. In comparison with the Fourier accelerated iter- 
ative schemes used for modeling lattice breakdown [23j . 
this algorithm significantly reduced the computational 
time required for solving large lattice systems. Using 
this numerical algorithm, we were able to investigate 
damage evolution in large (L = 512 for spring model) 
initially fully intact discrete lattice systems. However, 
due to insufficient number of available sample configura- 
tions, in this paper, we consider results up to L = 256 for 
spring models. For many lattice system sizes, the number 
of sample configurations, N con fi g , used are excessively 
large to reduce the statistical error in the numerical re- 
sults (see Table 1). In Table 1, the fraction of broken 
bonds (or damage density) for each of the lattice system 
sizes is obtained by dividing the number of broken bonds 
with the total number of bonds, N e i, present in the fully 
intact lattice system. For triangular lattice topology, 
N e i = (3L + 1)(L + 1). The lattice system sizes considered 
in this work are L = {8,16,24,32,64,128,256}. How- 
ever, since corrections to the scaling laws are strongest 
for small lattice systems, in the following, we use lattice 
sizes L > 16 for obtaining the scaling exponents. Table 
1 presents mean and standard deviations in the broken 
bond density (fraction of broken bonds) at the peak load 
and at failure for various triangular lattice system sizes. 
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FIG. 1: Envelope of a typical force-displacement response 
obtained using the RSM. 



III. DIFFUSIVE DAMAGE AND 
LOCALIZATION 

Qualitatively, damage evolution as described by break- 
ing of bonds is controlled by two competing aspects: dis- 
order and stress concentration in the vicinity of crack 
tips. In the case of strong disorder, bond breaking events 
occur in an uncorrelated manner in the initial stages of 
damage evolution and thus resemble percolation. As the 
damage starts to accumulate, some degree of correlation 
can be expected due to the presence of stress concentra- 
tion at the crack tips. A natural question to ask con- 
cerns the relevance of these correlations as failure is ap- 
proached. If correlations are irrelevant one should ob- 
serve percolation scaling up to failure, as in the case of 
infinite disorder. On the other hand, in the weak disor- 
der case, the current enhancement at the crack tips is so 
strong that a spanning crack is nucleated soon after a few 
bonds (or even a single bond) are broken Q. The inter- 
esting situation corresponds to the diffuse damage and 
localization regime, where a substantial amount of dam- 
age is accumulated prior to failure. Figure [21 presents 
the snapshots of damage evolution in a typical RSM sim- 
ulation of size L — 256. 

In order to investigate the localization of damage prior 
to failure, we divided the load-displacement response of a 
typical RSM simulation into 12 segments, with six equal 
segments each before and after the peak load. Figures 
13 and 01 present the snapshots of damage profiles within 
each segment of load-displacement curve of a typical sim- 
ulation with uniform threshold disorder for L = 256. 
Based on Figs. |21 and 01 it is clear that localization of 
damage occurs in the RSM prior to failure even for strong 
but finite disorder. In fact, the damage is diffusive in the 
initial stages of loading and extends upto almost the peak 
load. Around the peak load, the damage starts to local- 
ize and ultimately leads to failure, and hence the final 
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breakdown event is very different from the initial precur- 
sors upto the peak load. Similar behavior is observed in 
the random thresholds fuse model with both uniform and 
power law thresholds distributions [22|. 

In order to obtain a quantitative description of the 
damage localization process it is necessary to average 
the damage profiles over different realizations. Since the 
localization of damage can occur anywhere along the y 
direction of the lattice, a simple averaging of the dam- 
age profiles would yield a flat profile irrespective of the 
individual profile shapes in a single realization. In this 
study, we average the damage profiles by first shifting 
the damage profiles by the center of mass of the dam- 
age and then averaging. Alternatively, one could average 
the magnitude of the Fourier transforms of individual 
damage profiles thereby retaining the frequency content 
of damage profiles. The Fourier method eliminates any 
artificial biasing associated with the shifting of the indi- 
vidual profiles in the real-space |22|. 

Figure presents the average damage profiles for the 
damage accumulated up to the peak load by first shifting 
the damage profiles by the center of mass of the damage 
and then averaging over different samples. The results 
presented in Fig. [3] indicate that although the average 
damage profiles at smaller lattice system sizes are not 
completely flat, they flatten considerably as the lattice 
system size is increased. We tend thus to attribute the 
apparent profile to size effects. Indeed, for large system 
sizes (e.g. L — 128 and 256), the results clearly show that 
there is no localization at the peak load. Consequently, 
the localization of damage is mostly due to the damage 
accumulated between the peak load and failure, i.e., the 
final catastrophic breakdown event. Figure presents 
the data collapse of the average damage profiles for the 
damage accumulated between the peak load and failure 
using a power law scaling. A perfect collapse of the data 
is obtained using the form 

(Ap(y,L))/{Ap(0)) = f(\y-L/2\/0, (1) 

where the damage peak scales as (Ap(Q)) = L~ - 37 and 
the localization length scales as £ ~ L a , with a = 0.65 
(see Fig. |SJ| . The profile shapes decay exponentially at 
large system sizes. We have also tried a simple linear 
scaling of the form (Ap(y, £))/(Ap(0)) = f((y-L/2)/L), 
but the collapse of the data is not very good. The result 
for the fuse model is similar: the profile also displays 
exponential tails and the exponent is found to be a = 0.8. 

IV. SCALING OF DAMAGE DENSITY 

It has been noted in the previous section that the final 
breakdown event is very different from the initial pre- 
cursors. Thus, we consider the scaling of the number of 
broken bonds at the peak load, n p , that excludes the last 
catastrophic event. In Fig. [7|we plot Tip &-S cl function of 
the lattice size N e i . The data displays a reasonable power 
law behavior n p ~ N^, with b = 0.92. The exponent 



b = 0.92 is in close agreement with the value obtained 
for random thresholds fuse model using both triangular 
(b = 0.93) and diamond (6 = 0.91) lattice topologies 
|2^. The difference between the RSM and RFM models 
is marginal and may be attributed to the results obtained 
from the smaller lattice sizes, where corrections to the 
fractal scaling may exist. However, we have noticed some 
systematic deviations from the scaling form n p ~ by 
plotting rip/N^i vs N e i. Since the exponent b is close to 
one, the data could be equally well fit by a linear law 
times a logarithmic correction n p ~ N e i / \og(N e i) as sug- 
gested in Ref. |24| (see the inset of Fig. 0). Both these 
fits imply that in the limit of large lattices the fraction 
of broken bonds prior to fracture vanishes. 



V. SCALING OF DAMAGE DENSITY 
PROBABILITY DISTRIBUTIONS 

Since the final breakdown event is very different from 
the initial precursors up to the peak load, in this section, 
we present the scaling of the cumulative probability dis- 
tributions for the fraction of broken bonds at the peak 
load. The cumulative probability distribution for the 
damage density at the peak load is defined as the proba- 
bility Tlplpb, L) that a system of size L reaches peak load 
when the fraction of broken bonds equals pb = j^- , where 
n;, is the number of broken bonds. Figure |S] presents the 
cumulative probability distribution for the damage den- 
sity at the peak load in the random thresholds spring 
model for various system sizes L. By simply plotting 
the distribution in terms of p„ = p " p ' > = ^z£sl 

where fi n and a n denote the mean and standard devi- 
ation of the number of broken bonds at peak load, and 
p p and A p denote the mean and standard deviation of 
fraction of broken bonds at peak load (see Table 1), we 
obtain a very good collapse of the cumulative probabil- 
ity distribution of the damage density at the peak load. 
Fig. EH shows that Il p (p, L) may be expressed in a univer- 
sal scaling form such that n p (p, L) — H p (p p ) for different 
system sizes L. A similar collapse has been performed for 
the random thresholds fuse model in the Ref. 22]. The 
inset in Fig. [5] presents the comparison of the cumulative 
damage density probability distributions in the random 
thresholds spring and fuse models. The excellent collapse 
of the data in the inset of Fig. [5] suggests that the cu- 
mulative probability distribution for the damage density 
at the peak load, H p (pb,L) = H p (p p ), may be universal. 
Finally, the collapse of the data in Fig. ^| indicates that 
a Gaussian distribution adequately describes n p . 

In Ref. [23 , we have also checked that the distributions 
at failure in the random thresholds fuse model obey es- 
sentially the same laws, i.e., H(p) = Hf(pf) = H p (p p ), 
where Hf(pf) is the probability that a system of size L 
fails when the fraction of broken bonds equals pb, and 
pf is the corresponding reduced variable at failure. How- 
ever, in the RSM, although a reasonable collapse of the 
cumulative probability distribution of damage density at 
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FIG. 2: Snapshots of damage evolution in a typical simulation of size L — 256. Number of broken bonds at the peak load 
and at failure are 13864 and 16695, respectively. (l)-(9) represent the snapshots of damage after n b bonds are broken. (1) 
n b = 5000 (2) n b = 10000 (3) n b = 12000 (4) n b = 13000 (5) n b = 14000 (just after peak load) (6) n b = 15000 (7) n b = 15500 
(8) n b = 16000 (9) n b = 16500 (close to failure) 
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FIG. 3: Snapshots of pre-peak damage profiles of a typical 
RSM simulation with uniform threshold distribution on a tri- 
angular lattice of size L = 256. The damage is uniform in the 
pre-peak regime. 





FIG. 4: Snapshots of post-peak damage profiles of a typi- 
cal RSM simulation with uniform threshold distribution on 
a triangular lattice of size L = 256. The damage is clearly 
localized in the post-peak regime. 



failure can be obtained, the cumulative distributions of 
damage density at peak load and at failure appear to 
be different. In particular, the distribution, TLf(pf), at 
failure is not adequately described by a Gaussian distri- 
bution. The inadequacy of a Gaussian distribution in 
the post-peak regime may indicate the presence of a rel- 
atively stronger localization in the RSM compared with 
the RFM. 



FIG. 5: Average damage profiles at peak load obtained by first 
centering the data around the center of mass of the damage 
and then averging over different samples. For each of the 
samples, the damage profile is evaluated as p(y) = .^W. , 

where n b (y) denotes the number of broken bonds in the y th 
section. 



VI. UNIVERSALITY OF FRACTURE 
STRENGTH DISTRIBUTION 



In this study, we start the numerical simulation with 
a fully intact lattice system. The fracture strength of 
such a system is denned as the stress corresponding to 
the peak load of the lattice system response. 

Figure lXTT a) presents the fracture strength density dis- 
tributions for random thresholds spring model using the 
standard Lognormal variable, £, defined as £ 



Ln(o f )-ri 

where at refers to the fracture strength defined as the 
peak load divided by the system size L, and i] and ( 
refer to the mean and the standard deviation of the loga- 
rithm of <Tf . In order to verify the universality of fracture 
strength distribution, the fracture strength distributions 
from [l8| corresponding to random thresholds fuse model 
(RFM) using triangular lattice systems with uniform dis- 
order are presented in Fig. Illf a) along with those corre- 
sponding to random thresholds spring model. In particu- 
lar, Fig. II If a) shows the data for different lattice system 
sizes, L, corresponding to (a) triangular spring lattice, 
L = {8,16,24,32,64,128} and (b) triangular fuse lat- 
tices of sizes L = {4,8,16,24,32,64,128}. In all, there 
are 13 plots in Fig. Illf a). and the excellent collapse of 
the data for various spring and fuse lattices clearly in- 
dicates the universality of the fracture strength density 
distribution. The results presented in Fig. Illf a) are lim- 
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FIG. 6: Data collapse of the average profiles for the damage 
accumulated between peak load and failure using a power law 
scaling. We have considered the damage profiles for L — 
{16,24,32,64,128,256} system sizes. The average has been 
performed after shifting by the center of mass. The profiles 
show exponential tails. 
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FIG. 7: Scaling of number of broken bonds at peak load for 
triangular random thresholds spring lattices. The scaling ex- 
ponent b = 0.92 is very close to the exponent obtained for 
random thresholds fuse network using triangular (b = 0.93) 
and diamond (b = 0.91) lattices. The difference could be at- 
tributed to small size effects. The number of broken bonds 
at peak load can also be fit by a linear function times a loga- 
rithmic correction by plotting n p /N e i as a function of N e i in 
a log- linear plot (inset). 



ited only up to a system size of L = 128 due to the 
availability of fewer sample configurations for larger lat- 
tice systems. In order to attain a good collapse of the 
data for the density distributions, it is necessary to con- 
sider many sample configurations. On the other hand, 
good collapse of the data for the cumulative distributions 
can be achieved using fewer number of sample configu- 




0.25 



FIG. 8: The cumulative probability distribution for the frac- 
tion of broken bonds at the peak load for triangular spring 
lattices of different system sizes. 



rations. Figure. Illf b) presents the cumulative fracture 
strength versus the standard Lognormal variable, £, for 
random spring and fuse lattice networks for system sizes 
up to L = 512. In particular, in Fig. Illf b). we plot 
the numerical simulation results of RSM for system sizes 
L = {8, 16, 24, 32, 64, 128, 256} along with those of RFM 
for system sizes L = {4, 8, 16, 24, 32, 64, 128, 256, 512}. In 
all, there are about 16 curves (7 for triangular spring lat- 
tices and 9 for triangular fuse lattices) in Fig. Illf b). and 
the excellent collapse of the data suggests universality 
of fracture strength distribution. In Ref. [lfij . we have 
also presented the collapse of the fracture strength dis- 
tribution for different lattice topologies (such as trian- 
gular and diamond), which is consistent with the notion 
of universality of fracture strength distribution. That is, 
P(o < of — where P(o < of refers to the cumu- 

lative probability of fracture strength o < Of, ^ is a uni- 
versal function such that < \P < 1, and £ = Lnl ^' 7 ^' ri 
is the standard Lognormal variable. 

Figures IT2T a) and I12f b) present the modified Gumbel 
and Weibull fits for the fracture strength distribution of 
triangular spring lattice network using 



.4 



fc 



1 



'f . 



Ln c 



for the modified Gumbel distribution, and 



A = m Ln [ — 



Ln c 



(2) 



(3) 



for the Weibull distribution. In Eqs. |2Jand|3J fc, S, c and 

m are constants, and A is defined as 
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Ln(l-P(o f )) 
L 2 



(4) 



where P(of denotes the cumulative distribution. From 
these figures, it is clear that fracture strength data for 
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FIG. 9: The collapsed cumulative probability distribution for 
the fraction of broken bonds at the peak load in the ran- 
dom thresholds spring model (RSM) using triangular lattices 
of different system sizes (L = 8, 16, 24, 32, 64, 128, 256) with 
uniform disorder when plotted as a function of the reduced 
variable p p — (p — p p )/A p . In the inset, a comparison be- 
tween the cumulative probability distributions of the frac- 
tion of broken bonds at the peak load is presented for the 
RSM and RFM. For the RSM, triangular lattices of sizes 
[L = 8,16,24,32,64,128,256), and for the RFM, triangular 
lattices of sizes (L = 16, 24, 32, 64, 128, 256, 512) are plotted. 
In the RFM case, collapse of cumulative probability distribu- 
tions at the peak load for different lattice topologies (triangu- 
lar and diamond) and different disorder distributions (uniform 
and power law) is presented in Ref. |22T| . 



different lattice system sizes does not collapse on to a 
single straight line as it should, if the data were to follow 
Eq. J2J or © ■ This indicates that neither modified Gum- 
bel nor Weibull distributions may represent the fracture 
strengths distribution accurately for the RSM. In the Ref. 
|18|. similar conclusion has been drawn for the fracture 
strengths distribution of RFM. 

On the other hand, in Fig. I12f c) , we test the Lognor- 
mal description for fracture strengths by plotting the in- 
verse of the cumulative probability, &~ 1 (P(<jf)), against 
the standard Lognormal variable, £. In the above de- 
scription, $(•) denotes the standard normal probability 
function. In particular, in the Fig. I12f c) we present 
the Lognormal fit for the cumulative fracture strength 
distributions obtained for random thresholds spring and 
fuse models (i.e., for all of the 16 curves in Fig. Illf b)). 
Once again, this figure clearly indicates that the frac- 
ture strength distribution obtained for different lattice 
system sizes collapses onto a single curve, albeit minute 
deviation from straight line behavior is evident. We have 
also used the normal distribution to collapse the frac- 
ture strength data of triangular spring and fuse lattice 
systems. Although the data collapse is reasonable, it is 
not as good as that of Lognormal distribution based on 
Kolmogorov-Smirnov goodness-of-fit test. 
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FIG. 10: Normal distribution fit for the cumulative probabil- 
ity distributions of the fraction of broken bonds at the peak 
load for triangular spring lattices of different system sizes L 
= {8, 16, 24, 32, 64, 128, 256}. 



VII. SIZE EFFECTS IN THE MEAN 
FRACTURE STRENGTH 

The mean fracture strength data for various random 
threshold spring lattice system sizes is presented in Ta- 
ble 2. In Ref. for the RFM, we have suggested 
a scaling form F peak = C L a + C\ for the peak load, 
where Cq and C\ are constants. Correspondingly, the 



mean fracture strength defined as /j,/ 



j^- , is given by 
(if — Cq L a ~ L + We have used the same scaling law 
for the random thresholds spring model as well, and the 
result presented in Fig. El indicates that the exponent a 
is approximately equal to 0.97, which is once again in con- 
sistent with the a = 0.96 obtained for RFM using both 
triangular and diamond lattice topologies. The inset in 
Fig. El presents a power law fit (if ~ L~™ that is con- 
sistent with a Weibull distribution for fracture strengths. 
From the nonlinearity of the plots in the inset of Fig. 1131 
it is clear that the mean fracture strength does not fol- 
low a simple power law scaling that is consistent with a 
Weibull distribution for fracture strengths. We have also 
verified that the mean fracture strength does not follow 
a scaling law of form fij = A + Bi Ln l ^ na ^ lSi con " 
sistcnt with a modified Gumbel distribution for fracture 
strengths E E H3 • 

Since a very small negative exponent (a — 1) is equiv- 
alent to a logarithmic correction, i.e., for (1 — a) << 
1, L Q_1 ~ (log(L))~v, an alternative expression for 
the mean fracture strength may be obtained as [if — 

(LogL)'t' T' wnere M/ an d c are constants that are re- 
lated to the constants Co and C\ . This suggests that the 
mean fracture strength of the lattice system decreases 
very slowly with increasing lattice system size, and scales 
as (if rs ji^jyyp, with ip ks 0.15, for very large lattice sys- 
tems. 
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FIG. 11: Universality of fracture strength distribution in 
the random thresholds spring and fuse models, (a) Frac- 
ture strength density distributions for triangular spring lat- 
tices (L = {8, 16, 24, 32, 64, 128}) and triangular fuse lattices 
(L = {4, 8, 16, 24, 32, 64, 128}) with uniform disorder, (b) Cu- 
mulative fracture strength distribution for triangular spring 
lattices (L = {8, 16, 24, 32, 64, 128, 256}) and triangular fuse 
lattices (L = {4, 8, 16, 24, 32, 64, 128, 256, 512}) with uniform 
disorder. The collapse of the data in random spring and fuse 
models suggests universality of fracture strength distribution. 
In the RFM case, the universality of fracture strength distri- 
butions with respect to different lattice topologies is presented 
in Ref. [Ts||. 
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VIII. AVALANCHES 
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The avalanche size distribution, once the last event 
is excluded, is a power law followed by an exponential 
cutoff at large avalanche sizes (see Fig. I14|) ;25|. The 
cutoff size so is increasing with the lattice size, so that 
we can describe the distribution by a scaling form 

P(s,N) = s-Tg( S /N D ^), (5) 

where D represents the fractal dimension of the 
avalanches and N = (3L + 1)(L + 1) is the number of 



FIG. 12: Probability distribution fits for fracture strengths at 
the peak load in a triangular spring lattice network for dif- 
ferent lattice system sizes L = {8, 16, 24, 32, 64, 128, 256}. (a) 
Modified Gumbel distribution (top) . (b) Weibull distribution 
(middle), (c) Lognormal distribution fit for all the 16 curves 
(see Fig. Illf b'O (bottom). Since the data for different lattice 
system sizes does not collapse onto a single curve, Weibull 
distribution may not be an adequate fit for representing frac- 
ture strengths in the RSM. On the other hand, the collapse of 
the data in the reparametrized Lognormal distribution fit sug- 
gests that the Lognormal distribution describes the fracture 
strength distribution adequately. 
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FIG. 13: Proposed scaling law for the mean fracture strength 
(Fpeak = CoL a +C\). (1) Triangular spring network (symbol: 
-+-): a — 0.97. (2) Triangular fuse network (symbol: -*- 
): ol — 0.956; The corresponding Weibull fit for the mean 
fracture strength is shown in the inset. Nonlinearity of the 
plots in the inset suggests that mean fracture strength does 
not follow a power law scaling consistent with the Weibull 
distribution. 

bonds. Figure ED presents the data collapse of the dis- 
tribution of avalanche sizes using the exponents r = 2.5 
and D = 1.1. 
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FIG. 15: Data collapse of the avalanche size distributions 
excluding the hnal catastrophic event. The exponents used 
for the collapse are r = 2.5 (the reference line has this slope) 
and D — 1.1. The distributions have been logarithmically 
binned to reduce fluctuations. 



Fig. 14 of Ref. H^])- The cumulative distribution of 
last avalanche sizes for the RSM and RFM is presented 
in Figs. I17f a) and (b) respectively. While the distribu- 
tion is approximately Gaussian for RFM as shown by the 
data collapse (almost linear) in the inset of Fig. I17f b). 
there appears to be significant nonlinearity in the data 
collapse of the plots in the inset of Fig. lT7T a). This sug- 
gests that normal distribution may not be an adequate 
fit for representing the distribution of last avalanche size 
in the RSM model. We notice here that the post-peak 
regime is different in the two models because the RSM 
can fail because of a loss in rigidity. In general, the signif- 
icantly different nature of the last avalanche with respect 
to the precursors is revealed both by the distribution type 
(Gaussian or power law) and by its characteristic value, 
scaling as 2b ~ 1.36 or D = 1.1. This difference reflects 
the fact that the last avalanche is a catastrophic event 
corresponding to unstable crack growth, while precursors 
reflect metastable crack growth and the two processes are 
different. 



FIG. 14: The distribution of avalanche sizes (without the last 
catastrophic event) for triangular spring lattices of different 
sizes. 

So far we have considered avalanche statistics integrat- 
ing the distribution over all the values of the forces upto 
the peak load, but the avalanche signal is not station- 
ary: as the force increases so does the avalanche size. 
In particular, the last avalanche is much larger than the 
others. Its typical size grows as s m = (n/ — n p ) ~ N b , 
with b ~ 0.68 (see Fig. 1160 . which is once again in good 
agreement with the b ~ 0.7 value obtained for RFM (see 



IX. SUMMARY 

In this study we investigated the universality of ran- 
dom thresholds spring and fuse models using large scale 
numerical simulations and large number of sample con- 
figurations to reduce the statistical error in the numer- 
ical results. For both models, we considered triangu- 
lar lattice topology with uniform disorder and increased 
the load quasistatically. We performed several statistical 
measures characterizing the fracture process that can be 
summarized as follows: 
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FIG. 16: The mean avalanche size of the last catastrophic 
event (s m = (n/ — n p )) scales as a power law of N e i. Once 
again, the scaling exponent b = 0.68 for RSM is similar to the 
scaling exponent b ~ 0.7 obtained for RFM using triangular 
and diamond lattices (see Fig. 14 of Ref. |2^|). 

1 . Damage localization: the process of localization is 
similar in the RFM and RSM. Damage is accumu- 
lated in a uniform manner up to the peak load and 
then suddently localizes leading to complete fail- 
ure. This process is described by damage profiles, 
that are basically flat until peak load and show 
a peak, with exponential tails, in the post peak 
regime. The collapse of the damage profiles implies 
some small differences in the exponets for the two 
models. 

2. Damage density: The number of broken bonds at 
failure or at peak load scales with the lattice size. 
We have shown that the behavior in RFM and RSM 
is very similar, but in both cases it is not possibile 
to distinguish a power law from a linear behavior 
with a logarithmic correction. 

3. Damage distributions: The distribution of broken 
bonds at peak load follows the normal distribution 
for RFM and RSM. 

4. Fracture strength: The fracture strength distribu- 
tion is found to be lognormal for both models and 
the mean fracture strength scales logarithmically. 

5. Avalanches: The integrated avalanche distributions 
follow a power law in both models. The results of 
the RSM yields an exponent r = 5/2 that is very 
close to the exponent observed in global load shar- 
ing fiber bundle model (FBM), while larger devia- 
tions are found in the RFM (i.e. r = 2.7 13). 

Thus in conclusion, we can state that RFM and 
RSM are qualitatively very similar: distributions have 
the same forms, localization proceeds in the same way, 




-4 -2 2 4 6 8 
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FIG. 17: The collapsed cumulative distribution of last 
avalanche, (a) RSM (b) RFM. The insets in each of these 
figures show how well the data can be represented by a nor- 
mal distribution fit. The presence of significant nonlinearity 
of the data in these insets suggests that normal distribution 
may not be an adequate fit for representing the distribution 
of last avalanche size for the RSM model, whereas it may be 
an adequate fit for the RFM model. 



avalanches are similar. The only diffirences can be found 
in small quantitative deviations in exponents. We can 
not rule out that these are due to differences in the finite 
size behavior of the models and that at large scales the 
behavior is the same. In addition, the rigidity mechanism 
present in the RSM and not in the RSM could explain 
some deviations in the post-peak regime. 
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TABLE I: Mean and standard deviation of damage density 
at the peak load and failure in the random thresholds spring 
model using triangular lattice network with uniform disorder 
distribution. N co „fi g denotes the number of configurations 
used in averaging the results for each system size. p p and pf 
denote the mean fraction of broken bonds in a lattice system 
of size L at the peak load and at failure, respectively. Simi- 
larly, A p and A / denote the standard deviation of the fraction 
of broken bonds at the peak load and at failure respectively. 
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TABLE II: Peak load in the random thresholds spring model 
using triangular lattice network with uniform disorder distri- 
bution. N con f ig denotes the number of configurations used in 
averaging the results for each system size. 
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